September 1992 



HLRZ-92-69 



Chern-Simons term in the 
4-dimensional SU(2) Higgs Model 



F. Karsch 1 ' 2 , M. L. Laursen 1 '^ 
T. Neuhaus 2 and B. Plache 2 



1 HLRZ, c/o KFA Julich, P.O. Box 1913, D-5170 Jiilich, Germany 
2 Fak. f. Physik, Univ. Bielefeld, D-4800 Bielefeld, Germany 



Abstract 

Using Seibergs definition for the geometric charge in SU(2) lattice gauge theory, 
we have managed to apply it also to the Chern-Simons term. We checked the 
periodic structure and determined the Chern-Simons density on small lattices L 4 
and L 3 x 2, 4 with L = 4, 6, and 8 near the critical region in the SU(2) Higgs 
model. The data indicate that tunneling is increased at high temperature. 



a Talk presented by M. L. Laursen at the Symposium on Lattice Field Theory, Lattice 92, 
Amsterdam September 15-19, 1992. 
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1 Introduction 



Some years back 't Hooft found that the baryon number and the lepton number are 
not conserved in the electroweak theory JI|. While the B — L symmetry remains 
unbroken due to the anomaly cancellation, B + L is no longer conserved. This 
socalled baryon number violation is caused by the nontrivial topological winding of 
the SU(2) gauge fields. The anomaly of the fermionic current relates the winding of 
the gauge fields and changes the baryon number by an amount 

Nf /"*2 [ 

B(t 2 ) - B{t x ) = -/- 2 / d?xtr[F^ v ] (1.1) 

where Nf is the number of families of quarks and leptons. In the axial gauge Aq = 
we can relate the change in the baryon number to the change in the Chern-Simons 
number 

B(t 2 ) - B(h) = N f [Ncs(t 2 ) - Ncsih)] (1.2) 
where the Chern-Simons number Ncs is 

N cs = -g^ / d 3 x€ ijk tr[Ai{djA k + -A,A k )]. (1.3) 

At zero temperature such processes are exponentially suppressed as exp( — 2n/aw), 
a ~ 1/30. This is because any gauge field configuration which changes the winding 
number has an action at least that of the barrier height 27r/aw 

At high temperatures which prevail in the early universe such an exponential 
suppression is absent, since the system can pass over the barrier classically. The 
only suppression factor is the Boltzmann factor exp(—/3E) where E is the barrier 
height, and this factor is close to one ref. 0. Any baryon asymmetry generated 
at the GUT scale will get washed out as the universe approaches the electroweak 
phase transition from above. If one assumes that the transition is of first order, 
together with CP violating processes and thermal non-equilibrium (provided by the 
expansion of the universe) one can explain the baryon asymmetry in the universe. 

There exists semiclassical solutions of the Yang-Mills Higgs fields which are be- 
lieved to be important in this scenario. These solutions are known as sphalerons. 
They are static, have finite energy, but are unstable. The sphaleron sits so to speak 
on top of the barrier and it has a baryonic charge of 1/2. While the instanton tunnels 
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between two Chern-Simons vacua (in the axial gauge) one must imagine a different 
time dependent solution which interpolates between a Chern-Simons vacuum and 
the top of the barrier. 

There are some lattice studies of baryon number violating processes in the 2d 
- Abelian Higgs model [Q], and in the 4d - SU(2) Higgs model [|J. The configu- 
rations are prepared at high temperature and the system is allowed to change via 
the classical Hamiltonian equation of motion. Since the axial gauge is used, the 
Gauss constraint must be implemented in addition. The Chern-Simons term Ncs{t) 
is monitored during the time evolution as a function of the temperature. When 
the system passes through a sphaleron transition one finds ANcs = ±1. All these 
calculations are done in the real time formalism. 

We have initiated work on the 4d - SU(2) Higgs model in Euclidean time, trying 
to see how the temperature influences tunneling. To evaluate N C s we have used 
a geometric definition given by Seiberg see ref. ||. It requires a two dimensional 
numerical integration, but we have an efficient vectorized code which is an improve- 
ment of the code used in ref. ||. See also ref. for alternative definitions of the 
Chern-Simons term. 

2 Topological charge and the Chern-Simons term 
in the continuum 

We will first define the topological charge in SU(2). The gauge field is and 
the gauge field tensor is = d^A u — d u A^ + L4 M ,A„]. Under a a local gauge 
transformation g the gauge field changes as: 5A^ = g'^A^ + d tl ]g(x), while the 
gauge field tensor transforms gauge covariantly F^ u — > g^F^g. The topological 
charge Q is gauge invariant and an integer, 



The manifold is denoted M and we shall assume that its boundary dM is a three 
sphere S 3 . The topological charge density q can be written as a perfect derivative 




(2.1) 
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(2.2) 



Q = ~ 
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where the Chern-Simons density is 

K » = -^^MMdpAa + \A p A a )\. (2.3) 

It is gauge variant and changes under the gauge transformation g by an amount 
{Q v = d v gg~ l ) 

8K» = - ^^e pupa tr[g u g p g a ] 

d u tr[G P A a ]. (2.4) 
We define the (timelike) Chern-Simons number Ncs as follows: 

N C s = I d 3 xK £ Z. (2.5) 

J dM 

While N C s is om y an integer for pure gauge configurations, the gauge variation is 
an integer (the boundary term vanishes) 



5N t 



cs 



^e 0vpa [ d 3 xtr[g u g p g a ] e z. (2.6) 

All 2 JdM 



This follows also from homotopy theory using the mapping g : S* 3 — > SU{2) = S 3 . 
Such mappings are characterized with the homotopy class n 3 (S' 3 ) G Z. 



3 Topological charge and the Chern-Simons term 
on the lattice 



We will now consider the lattice version of the topological charge and the Chern- 
Simons number. We will use a geometric definition. Problems with dislocations will 
be ignored here. The manifold is a four torus M = T 4 and we will cover M with 
cells (hypercubes) c(n). Let the gauge potential A" be defined on c(n) and likewise 
A^i 1 be defined on c(n — ft). At the faces (cubes) f(n, fx) = c(n — p) fl c(n), we can 
relate the two potentials by a transition function v U)il 

= <»\K + d u ]v ntil . (3.1) 

In Liischers version one first fixes to a local complete axial gauge in each c(n). This 
will define v n>fl at the corners of the hypercube. It is then possible to extend it to 
the whole cube. The topological charge is (S v = s^d^s, V v = p^d^p), 

Q L = E9 L (») = H(- l T(K, - k n+ p,p), 

n n,n 
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+ f d 2 xtr[V p S a ]. (3.2) 

87T Jdf 

The function s is denned on the cube, while p is defined on the boundary of the cube. 
The actual expressions are given in ref. |J. In Seibergs version no local gauge fixing is 
performed, but otherwise the same interpolation is performed. Replace (s,p, k njfl ) — > 
(S, P, K n ^). The difference is that that S and P only depend on the original gauge 
fields in the cube. Then N C s = J2n s Kn s ,n is nothing but a Chern-Simons term (the 
summation is over the spatial lattice only). Like in the continuum it is only an integer 
for pure gauge field configurations, but under gauge transformations it changes by an 
integer. The corresponding topological charge is defined as Q s = J2nQ S { n )y where 
— 1/2 < q s (n) < 1/2. By restricting the charge to this interval we will have a gauge 
invariant charge definition. We find that both the topological charge and the Chern- 
Simons term have the correct naive continuum limit. As an interesting corollary we 
verify that the two charge definitions are related. If we introduce the section of the 
Liischer bundle: w(x), x G dc(n), then s(x) = w(0)S(x)w^ 1 (x) Inserting this in the 
above expressions yields q s (n) = q L (n)—q w (n) where q w (n) is the topological charge 
(integer) of the section. Therefore q s (n) = q L (n) up to integers ||. For smooth fields 
like instantons they always agree, while for realistic configurations this is true for 
almost every hypercube. It is interesting to look at the Chern-Simons number at 
each time slice of an instanton in the complete axial gauge. This is shown in Fig. |l|. 
We start in the vacuum sector Ncs = and move towards the other vacuum sector 
Ncs = 1- On the last time slice we must return to the first vacuum due to periodic 
boundary conditions. 



4 Tests and Monte-Carlo results for the Chern- 
Simons density 

We have used an SU(2) Higgs model with action: 

S= ~ § E tr[U n ^] + \J2(^n - 1) 2 

n,fi<v n 

- Kj2tr[& n U n ,^ n+tl }. (4.1) 
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We always tried to work close to the Higgs phase transition. We looked at the 
periodic structure for a 4 3 x 2 lattice at (f3, k, A) = (2.25,0.3,0.5). We did 1000 
configurations and we measured Ncs without any gauge fixing. In Fig. |2| we have 
plotted the Chern-Simons probability. The periodic structure is obvious. Notice 
that most of the configurations have Ncs close to an integer, and therefore can 
be interpreted as being pure gauge. We checked that Ncs indeed changed by an 
integer under an axial gauge transformation. Next we compared the probabilities 
for the two lattices 6 4 and 6 3 x 2. For these lattices we used 50 Landau gauge fixing 
sweeps, to make the integrals converge fast. In both cases we have 6000 Chern- 
Simons numbers. The results for the restricted and therefore gauge invariant Ncs 
( — 1/2 < Ncs < 1/2) are shown in Fig. ^| and Fig. |j. There is a trend in the 
direction of a flatter distribution at finite temperature. We take this as evidence 
that the system tunnels more often. A few words about the integrations done. One 
of the integrations can be done analytically so we are left with a two dimensional 
integral. We have used the following strategy, which turned out to be quite efficient. 
Perform a Gaussian integration with 8x8 points and store the results for the eight 
Kn^'s in each hypercube. Redo the same thing with 16 x 16 points and compare 
the results for each K n ^. If the relative difference is less than 0.0001 we accept 
the contribution. Otherwise we collect the K n ^s which have not yet converged. 
These we redo with 32 x 32 points instead. Compare with the previous values and 
repeat the procedure with 64 x 64 points. Usually, at this point only a few K^'s 
have not converged, so for these we use a library integration routine with interval 
adaption. The typical time for one topological charge on a 6 4 lattice is 100 seconds 
on the CRAY-YMP. The charges are integers up to errors of the order 10 -4 . The 
corresponding time for one Chern-Simons number in a timeslice is around 5 seconds. 
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Figure 1: Profile of Ncs through an instanton configuration. 
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Figure 2: The periodic structure for Ncs- 
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Figure 3: The Chern-Simons probability at zero temperature as a function of Ncs- 



Figure 4: The Chern-Simons probability at finite temperature as a function of Ncs- 
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